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Abstract. This is the first attempt to model the kinematics of 
a CME launch and the resulting EUV dimming quantitatively 
with a self-consistent model. Our 4D-model assumes self- 
similar expansion of a spherical CME geometry that consists 
of a CME front with density compression and a cavity with 
density rarefaction, satisfying mass conservation of the to- 
tal CME and swept-up corona. The model contains 14 free 
parameters and is fitted to the 2008 March 25 CME event 
observed with STEREO/A and B. Our model is able to re- 
produce the observed CME expansion and related EUV dim- 
ming during the initial phase from 18:30 UT to 19:00 UT. 
The CME kinematics can be characterized by a constant ac- 
celeration (i.e., a constant magnetic driving force). While the 
observations of EUVI/A are consistent with a spherical bub- 
ble geometry, we detect significant asymmetries and density 
inhomogeneities with EUVI/B. This new forward-modeling 
method demonstrates how the observed EUV dimming can 
be used to model physical parameters of the CME source re- 
gion, the CME geometry, and CME kinematics. 
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1 Introduction 

It has become a generally agreed concept that the EUV dim- 
ming observed during the onset of a coronal mass ejection 
(CME) manifests the coronal mass loss of the CME, and 
thus we basically expect a one-to-one correlation between 
the detections of CMEs and EUV dimmings, unless there ex- 
ist special circumstances. For instance, the CME could orig- 
inate behind the limb, in which case the EUV dimming is 
obscured, or the CME could start in the upper corona, where 
there is little EUV emission because of the gravitational strat- 
ification. The latter case would imply very low masses com- 
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pared with a CME that originates at the base of the corona, 
i.e., « 10% at two thermal scale heights. However, there 
exists a case with an average CME mass that did not leave 
any footprints behind in EUV (Robbrecht et al. 2009). A 
statistical study on the simultaneous detection of EUV dim- 
mings and CMEs has recently been performed by Bewsher 
et al. (2008). This study based on SOHO/CDS and LASCO 
data confirms a 55% association rate of dimming events with 
CMEs, and vice versa a 84% association rate of CMEs with 
dimming events. Some of the non-associated events may be 
subject to occultation, CME detection sensitivity, or incom- 
plete temperature coverage in EUV and soft X-rays. Perhaps 
the CME-dimming association rate will reach 100% once the 
STEREO spacecraft arrive at a separation of 180° and cover 
all equatorial latitudes of the Sun. 

A number of studies have been carried out by using the de- 
tection of coronal dimming to identify CME source regions, 
focusing on transient coronal holes caused by filament erup- 
tions (Rust 1983; Watanabe et al. 1992), EUV dimming at 
CME onsets (Harrison 1997; Aschwanden et al. 1999), soft 
X-ray dimming after CMEs (Sterling & Hudson 1997), soft 
X-ray dimming after a prominence eruption (Gopalswamy 
& Hanaoka 1998), simultaneous dimming in soft X-rays and 
EUV during CME launches (Zarro et al. 1999; Harrison & 
Lyons 2000; Harrison et al. 2003), determinations of CME 
masses from EUV dimming from spectroscopic data (Har- 
rison & Lyons 2000; Harrison et al. 2003) or from EUV 
imaging data (Zhukov and Auchere 2004; Aschwanden et 
al. 2009b). All these studies support the conclusion that dim- 
mings in the corona (either detected in EUV, soft X-rays, or 
both) are unmistakable signatures of CME launches, and thus 
can be used vice versa to identify the mutual phenomena. 

In this study here we attempt for the first time to model 
the kinematics of a CME and the resulting EUV dimming 
quantitatively, which provides us unique physical parameters 
of the CME source region and on the CME kinematics in the 
initial acceleration phase. 
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2 Model Assumptions 

The dynamics of a CME can often be characterized by a rapid 
expansion of a magnetically unstable coronal volume that ex- 
pands from the lower corona upward into the heliosphere. 
Different shapes have been used to approximately describe 
the 3D geometry of a CME, such as a spherical bubble, an 
ice-cone, a crescent, or a helical flux rope, which expand in 
a self- similar fashion and approximately maintain the aspect 
ratio in vertical and horizontal directions during the initial 
phase of the expansion. Here we develop a four-dimensional 
(4D=3D+T) model that describes the 3D evolution of the 
CME geometry in time (T) in terms of 4D electron density 
distributions n e (x,y, z,t) that allow us also to predict and 
forward-fit a corresponding EUV intensity image data cube 
Ieuv(%, t) in an observed wavelength. 

For the sake of simplicity we start in our model here with 
the simplest case, assuming: (1) spherical 3D geometry for 
the CME front and cavity; (2) self- similar expansion in time; 
(3) density compression in the CME front and adiabatic vol- 
ume expansion in the CME cavity; (4) mass conservation 
for the sum of the CME front, cavity, and external coronal 
volume; (5) hydrostatic (gravitational stratification) or super- 
hydrostatic density scale heights; (6) line-tying condition for 
the magnetic field at the CME base; and (7) a magnetic driv- 
ing force that is constant during the time interval of the initial 
expansion phase. This scenario is consistent with the tradi- 
tional characterization of a typical CME morphology in three 
parts, including a CME front (leading edge), a cavity, and a 
filament (although we do not model the filament part). The 
expanding CME bubble sweeps up the coronal plasma that 
appears as a bright rim at the observed "CME front" or lead- 
ing edge. The interior of the CME bubble exhibits a rapid 
decrease in electron density due to the adiabatic expansion, 
which renders the inside of the CME bubble darker in EUV 
and appears as the observed "CME cavity". The assumption 
of adiabatic expansion implies no mass and energy exchange 
across the outer CME boundary, and thus is consistent with 
the assumption of a low plasma ^-parameter in the corona 
with perfect magnetic confinement, while the CME bubble 
will become leaking in the outer corona and heliosphere, 
where the plasma ^-parameter exceeds unity (not included 
in our model here). 



3 Analytical Model 

A spherical 3D geometry can be characterized by one single 
free parameter, the radius R of the sphere. The self-similar 
expansion maintains the spherical shape, so the boundary of 
the CME bubble can still be parameterized by a single time- 
dependent radius R(t). The time-dependence of the CME 
expansion is controlled by magnetic forces, e.g., by a Lorentz 
force or hoop force. For sake of simplicity we assume a con- 
stant force during the initial phase of the CME expansion, 
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Fig. 1. The definition of the self-similar geometry of the CME 
model is depicted for four different times, consisting of a cyclidrical 
base and a spherical shell. The height of the centroid of the spher- 
ical CME volume is h(t), the outer radius of the CME sphere is 
R(t), and the inner radius of the CME front is r(t). These param- 
eters increase quadratically with time. The circular footpoint area 
of the CME with radius ro stays invariant during the self- similar 
expansion in order to satisfy the line-tying condition of the coronal 
magnetic field at the footpoints. 

which corresponds to a constant acceleration clr and requires 
three free parameters (i? , vr, clr) to characterize the radial 
CME expansion, 

R(t)=Ro + v R (t-to) + ^a R (t-t ) 2 for*>t , (1) 

where Ro = R(t = to) is the initial radius at starting time 
to, vr is the initial velocity and clr is the acceleration of the 
radial expansion. 

For the motion of the CME centroid at height h(t) we as- 
sume a similar quadratic parameterization, 

h(t) = h + v h (t - to) + ^a h (t - to) 2 for t > t , (2) 

where ho = h(t = to) is the initial height at starting time 
to, Vh is the initial velocity and ah is the acceleration of the 
vertical motion. This parameterization is consistent with a fit 
to a theoretical MHD simulation of a breakout CME (Lynch 
et al. 2004) as well as with kinematic fits to observed CMEs 
(Byrne et al. 2009). 

Further we constrain the CME geometry with a cylindri- 
cal footpoint area of radius ro, which connects from the so- 
lar surface to the lowest part of the spherical CME bubble. 
In order to ensure magnetic line-tying at the footpoints, the 
CME bubble should always be located above the cyclidrical 
footpoint base, which requires that the initial height satisfies 
ho > ro and the acceleration constants are ah > a r . We 
visualize the model geometry in Fig. 1. 

The time-invariant CME footprint area allows us a simple 
mass estimate of the CME from the cylindrical volume inte- 
grated over a vertical scale height, since the spherical CME 
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bubble will eventually move to large heights with no addi- 
tional mass gain (at time t n > to, see right-hand panel in 
Fig.l). 

Assuming adiabatic expansion inside the CME cavity, the 
electron density in the confined plasma decreases recipro- 
cally to the expanding volume, i.e., 

n e (t) = n (^j 3 , (3) 

so it drops with the third power as a function of time from 
the initial value no (of the average density inside the CME). 

For the mass distribution inside the CME we distinguish 
between a compression region at the outer envelope, con- 
taining the CME front, and a rarefaction region in the in- 
side, which is also called CME cavity. We define an aver- 
age width wo of the CME front that is assumed to be ap- 
proximately constant during the self- similar expansion of the 
CME. While the radius R(t) demarcates the outer radius of 
the CME front, we denote the inner radius of the CME front 
or the radius of the cavity with r(t), 



r(t) = R(t) - wo , 



(4) 



which has an initial value of ro = r(t = to) = Ro — wo. 
The total volume Vo of the CME is composed of a spherical 
volume with radius R(t) and the cylindrical volume beneath 
the CME with a vertical height of (ho — ro), 



V(t) = -7TR(t) 3 +7TR 2 [h(t)-R(t)], 

which has an initial volume value of Vo, 

Vo = V(t = t ) = ^nR% + nRl[h - R ] ■ 

The volume of the CME front Vp{t) is essentially the dif- 
ference between the outer and inner sphere (neglecting the 
cylindrical base at the footpoint) 



(5) 



(6) 



V F (t) « ^R(tf - ^r(i) 3 , 
while the volume Vc of the cavity is, 



(7) 



V c (t) « -7rr(i) 3 . (8) 

We have now to define the time-dependent densities in the 
CME, for both the CME front, which sweeps up plasma dur- 
ing its expansion, as well as for the CME cavity, which rar- 
ifies due to the adiabatic expansion. The total mass m E (i) 
of the plasma that is swept up from the external corona in a 
CME corresponds to the total CME volume V(t) minus the 
initial volume of the CME cavity, 



m E (t) =m p <n E > [V(t) -V c (t = t )] 
= m p < n E > ^7r[R(t) 3 - r§] 



(9) 



where m p is the mass of the hydrogen atom and < n E > is 
the electron density in the external corona averaged over the 



CME volume. The same mass has to be contained inside the 
volume Vf of the CME front, 



m E (t) = m p < tif > Vp(t) 

= m p < n F > §7r[i?(t) 3 



r(tf 



(10) 



Thus, mass conservation yields a ratio of the average electron 
density < n E > in the CME front and the average external 
density < n E > of 



Qn, front (t) — 



< n F > _ R(t) 3 -rjj 



< n E > R(t) 3 -r(t) 3 



(11) 



This density ratio amounts to unity at the initial time, i.e., 
Qn(t = to) = 1 and monotonically increases with time. The 
maximum value of the density jump in MHD shocks derived 
from the Rankine-Hugoniot relations (e.g., Priest 1982) is 
theoretically q n ,raax = 4. Fast CMEs are expected to be 
supersonic and will have a higher compression factor at the 
CME front than slower CMEs. Thus we keep the maximum 
compression factor q n , m ax as a free parameter, keeping in 
mind that physically meaningful solutions should be in the 
range of 1 < q n , m ax ~ 4. Since we prescribe both the width 
wo of the CME front as well as a maximum density compres- 
sion ratio q n ,max we obtain a definition of the critical value 
p(t) for the cavity radius r(t) when the prescribed maximum 
density compression q n ,max is reached (using Eq. 11), 



q n 



Ritf-rl 



R{tf- P {tf ' 
which yields the critical radius p(t), 



p(t) 



R(t) 



3 mf-ri) 

qn,max 



1/3 



(12) 



(13) 



Therefore, only plasma outside this critical radius p(t) can 
be compressed in the CME front, while the plasma inside 
this critical radius dilutes by adiabatic expansion and forms 
the cavity, yielding an average density ratio q n ,cav inside the 
cavity (according to Eq. 3), 



qn,cav (t) 



M _ f[ro/p(t)} 3 for p(t)>r(t) 



no * \[ro/r(t)} 3 for p(t) < r(t) 



.(14) 



Our numerical model of a spherical CME expansion has a 
total of 14 free parameters: 3 positional parameters (the he- 
liographic coordinates (/, b) and height ho of the initial CME 
centroid position), 5 kinematic parameters (starting time to, 
velocities Vh, vr, accelerations a^, or), 2 geometric param- 
eters (initial radius ro and width wo of the CME front), and 
4 physical parameters (coronal base density no, maximum 
density compression factor q n ,max in the CME front, the 
mean coronal temperature To at the observed wavelength 
filter), and a vertical density scale height factor (or super- 
hydro staticity factor) q\ that expresses the ratio of the ef- 
fective density scale height to the hydrostatic scale height 
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Fig. 2. Partial EUVI base-difference images from STEREO/A 171 A during the time interval of 2008 Mar 25, 18:36-18:56 UT (with the 
initial image at 18:36 UT subtracted). The envelope of the lateral CME front is fitted with a circle. The dashed lines mark the locations 
where cross-sectional brightness profiles are extracted for fitting, shown in Fig. 5. 



at temperature To). The temperature To defines the hydro- 
static scale height At of the corona in the observed tempera- 
ture range, which enters the definition of the effective density 
scale height A (e.g., Eq. [3.1.16] in Aschwanden 2004), 



A = A T qx = 47 



T f 



o 



1 MK 



q x [Mm] 



(15) 



Thus, assuming an exponentially stratified atmosphere 
(Eq. 15), a density compression factor q n (t) < q n ,max in 



the CME front (Eq. 12), and adiabatic expansion inside the 
CME cavity (Eq. 14), we have the following time-dependent 
3D density model: 



n e (h, t) = no exp (— h/X)x 

'1 for d > R(t) 

x < q n jront(t) for r(t) < d < R(t) 
(t) for d < r(t) 



(16) 
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18:36:00 




Fig. 3. EUVI base-difference images from STEREO/B 171 A in similar representation as for STEREO/A in Fig. 2. Note that the CME 
origin is in the front of the solar disk in EUVI/B, while it is behind the east limb for EUVI/A. 



where d is the distance of an arbitrary location with 3D 
coordinates (x,y,z) to the instantaneous center position 

[xo(t), yo(t), zo(t)] of the CME, 



4 Forward-Fitting of Model to Observations 
4.1 STEREO/EUVI Observations 



d=y/[x- x (t)]i + [y- yo (t)] 2 + [z ~ z (t)] 2 , (17) 



which is located at height h(t) vertically above the helio- 
graphic position (Z, b). 



One CME event observed with STEREO that appears as a 
spherically expanding shell, and thus is most suitable for fit- 
ting with our analytical model, is the 2008-Mar-25, 18:30 
UT, event. This CME occurred near the East limb for space- 
craft STEREO/Ahead, and was observed on the frontside 
of the solar disk with spacecraft STEREO/Behind. Some 
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preliminary analysis of this event regarding CME expan- 
sion and EUV dimming can be found in Aschwanden et al. 
(2009a), the CME mass was determined in white light with 
STEREO/COR-2 (Colaninno and Vourlidas 2009) and with 
STEREO/EUVI (Aschwanden et al. 2009b), and the 3D ge- 
ometry was modeled with forward- fitting of various geomet- 
ric models to the white-light observations (Thernisien et al. 
2009; Temmer et al. 2009; Maloney et al. 2009; Mierla et al. 
2009a,b). While most previous studies model the white-light 
emission of this CME, typically a few solar radii away from 
the Sun, our model applies directly to the CME source re- 
gion in the lower corona, as observed in EUV. We follow the 
method outlined in Aschwanden et al. (2009a). Our model 
also quantifies the geometry and kinematics of the CME, as 
well as the EUV dimming associated with the launch of the 
CME. 

4.2 Fitting of Positional Parameters 

Figures 2 and 3 show sequences of 16 (partial) EUV images, 
simultaneously observed with STEREO/A and B with a ca- 
dence of 75 s during the time interval of 18:3 6- 18:56 UT on 
2008-Mar-25. In order to determine the positional param- 
eters of the CME as a function of time we trace the outer 
envelope of the CME bubble (by visual clicking of 3 points) 
in each image and each spacecraft and fit a circle through 
the 3 points in each image. The selected points for fitting 
the position of the CME bubble were generally chosen in the 
brightest features of the lateral CME flanks, but could not al- 
ways been traced unambiguously in cases with multiple flank 
features. In those cases we traced the features that were clos- 
est to a continuously expanding solution. The radii and y- 
positions of the circular fits are fully constrained from the 
STEREO/A images, so that only the x-positions of the cen- 
troid of the spherical shell need to be fitted in the epipolar 
STEREO/B images. We note that the fits of the CME bubble 
roughly agree with the envelopes of the difference flux in the 
STEREO/B images initially (up to 18:48 UT), while there is 
a discrepancy later on. Apparently the CME has a more com- 
plex geometry than our spherical bubble model, which needs 
to be investigated further. 

This procedure yields the CME centroid positions 
[xA(ti),y A (U)} and [x B (ti),y B (U)] for the time sequence 
U,i = 1, 16. The images in Fig. 2 and 3 are displayed 
as a baseline difference (by subtracting a pre-CME image at 
18:36 UT) to enhance the contrast of the CME edge. The 
circular fits to the CME outer boundaries are overlayed in 
Fig. 2 and 3. Both images have been coaligned and ro- 
tated into epipolar coordinates (Inhester et al. 2006), so that 
the y-coordinates of a corresponding feature are identical in 
the spacecraft A and B images, while the x-coordinates dif- 
fer according to the spacecraft separation angle a sep , which 
amounts to a sep = 47.17° at the time of the CME event. 
The epipolar coordinates measured from both spacecraft are 



then related to the heliographic longitude I, latitude b, and 
distance r c from Sun center as follows, 

xa = t c cos(b) sin(Z^) 

Va = Vb = r c sin(6) , (18) 

x B = r c cos(b) sm(l A + a sep ) 

which can directly be solved to obtain the spherical (epipo- 
lar) coordinates (I a, 6, r c ), 

I A = a rctan (sin (a sep ) / \x A /x B - cos (a 3ep )) 

r c = y/y 2 A + [x A /sm(l A )]* (19) 

b = arcsin (yA/r c ) 

Therefore, using stereoscopic triangulation, we can directly 
determine the spherical coordinates (Z^, bi, r c ^), i = 1, 16 
for all 16 time frames, as well as obtain the CME curva- 
ture radii R(U) from the circular fits to the CMEs. We plot 
the so obtained observables lA(i), b(t), R(t), and h(t) = 
r c (t) — Rq in Fig. 3 and determine our model parame- 
ters I and b from the averages. We obtain a longitude of 
I A = -102.4° ± 0.9° (for spacecraft STEREO/A), l B = 
I a + a S e P = -54.9° ± 0.9° (for spacecraft STEREO/B), 
and a latitude b = -8.8° ± 0.6°. Thus, the CME source 
region is clearly occulted for STEREO/A. These epipolar 
coordinates can be rotated into an ecliptic coordinate sys- 
tem by the tilt angle /3ab = 2.66° of the spacecraft A/B 
plane. Viewed from Earth, the longitude is approximately 
l E « -102.4 + a sep /2 « -78.8°. Thus, the CME source 
region is 12° behind the East limb when seen from Earth. 
This explains why the EUV dimming is seen uncontaminted 
from post-flare loops, which are seen by STEREO/B but hid- 
den for STEREO/A. 

4.3 Fitting of Kinematic Parameters 

We plot also the observables h(t) and R(t) in Fig. 4 and 
determine the model parameters ho, Ro, ah, clr by fitting the 
quadratic functions R(t) (Eq. 1) and h(t) (Eq. 2), for which 
we obtain the starting time t = 18.64 hrs (18:38 UT), the 
initial CME radius Ro = 45 Mm, the initial height ho = 
45 Mm, and the accelerations cir = 0.54 km s -2 for the 
CME radius expansion, and ah = 0.52 km s -2 for the height 
motion. The initial velocity is found to be negligible (vr « 
and Vh ~ 0). We estimate the accuracy of the acceleration 
values to be of order w 10%, based on the uncertainty of 
defining the leading edge of the CME. Thus, we determined 
9 out of the 14 free parameters of our model sofar. 

Note that the acceleration measured here refers to the 
very origin of the CME in low coronal heights of < 0.6 
solar radii observed in EUVI data. The acceleration is ex- 
pected to be initially high and to rapidly decline further 
out, when the driving magnetic forces decrease at large al- 
titudes. This explains why our values for the acceleration 
in low coronal heights are significantly higher than mea- 
sured further out in the heliosphere, typically in the or- 
der of tens of m s _1 in height ranges of 5-22 solar radii, 
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Fig. 4. Measurements of the heliographic longitude Za(£) (top 
panel) and latitude b(t) (second panel) of the CME centroid, the 
radius R(t) of the CME sphere (third panel), and CME centroid 
height h(t) (bottom). The average values are I a — —102.4° ± 0.9° 
and b = -8.8° ± 0.9° for the heliographic position. The CME ra- 
dius R(t) and height h(t) are fitted with quadratic functions (dashed 
curves), yielding the constants t = 18.64 hrs (18:38 UT), R = 45 
Mm, ho = 45 Mm, and accelerations a# = 0.54 km s -2 and 
ah = 0.52 km s -2 (see Eqs. 1-2). 



as measured with SOHO/LASCO. SOHO/LASCO reported 
even a slightly negative acceleration at altitutes of 5-22 solar 
radii. The driving magnetic forces that accelerate a CME are 
clearly confined to much lower altitudes. 

4.4 Fitting of Geometric Parameters 

We model the 3D geometry of the CME bubble with the 
time-dependent radius R(t) and the width w of the CME 
compression region. In Fig. 5 we show cross-sectional 
EUV brightness profiles across the CME in horizontal di- 
rection (parallel to the solar surface) and in vertical direction 
for the EUVI/A 171 A observations (indicated with dotted 
lines in Fig. 2). These baseline- subtracted profiles clearly 
show a progressive dimming with a propagating bright rim at 
the CME boundary, which corresponds to the density com- 
pression region at the lateral expansion fronts of the CME. 
The bright rims are clearly visible in the images during 



18:46—18:56 UT shown in Fig. 2. The average width of the 
observed bright rims is wq « 10 Mm, a value we adopt in 
our model. 

4.5 Fitting of Physical Parameters 

Finally we are left with the four physical parameters 
To, q\, no, and q n ,max- Since we show here only data ob- 
tained with the 171 A filter, the mean temperature is con- 
strained by the peak temperature of the instrumental EUVI 
response function, which is at To = 0.96 MK. This con- 
strains the thermal scale height to At = 47, 000 x 0.96 = 
45,000 km. 

The remaining free parameters q\, an d q n ,max need 
to be determined from best-fit solutions by forward- fitting of 
simulated EUV brightness images (or profiles, as shown in 
Fig. 5) to observed EUV brightness images (or profiles). The 
EUV emission measure in each pixel position (x, y) can be 
computed by line-of-sight integration along the z-axis in our 
3D density cube n e (x, y, z) per pixel area dA for each time 



EM(x,y) = J n 2 e (x,y,z)dzdA , 



(20) 



from which the intensity Im{x,y) in the model image in 
units of DN s _1 can be directly obtained by multiplying with 
the instrumental response function R 171 (T) of the 171 A fil- 
ter, 



hri(x,y) = J EM(x,y)R 171 (T)dT 
^EM(x,y)R 171 (T )AT 171 



(21) 



where R 171 = 3632 x 10~ 44 DN s" 1 cm 3 MK" 1 and the 
FWHM of the 171 filter is AT m = 1.25 - 0.51 = 0.74 
MK. 

In Fig. 6 we show best-fit solutions of horizontal and ver- 
tical brightness profiles. The absolute flux level is propor- 
tional to the coronal base density squared, which we obtain 
by minimizing the mean flux difference between simulated 
and observed flux profiles. We obtain a best-fit value of 
no = 6.5 x 10 8 cm -3 . The super-hydro staticity factor is 
most sensitive to the vertical flux profile (Fig. 5, right-hand 
side panels), for which we find a best-fit value of #a = 1-45. 
Thus, the average density scale height in the CME region 
is slightly super-hydrostatic, as expected for dynamic pro- 
cesses. These values are typical for quiet-Sun corona and 
active region conditions (see Figs. 6 and 10 in Aschwanden 
and Acton 2001). 

The last free parameter, the maximum density compres- 
sion factor q n ,max, affects mostly the brightness of the CME 
rims. Fitting the brightness excess at the CME rims at those 
times where bright rims are visible are consistent with a value 

of qn,max ~ 2. 
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Fig. 5. Cross-sectional (baseline- subtracted) EUV brightness profiles in horizontal direction above the limb (left-hand panels) and in vertical 
direction through the centroid of the dimming region (right-hand panels), as marked with dashed lines in Fig. 2, for the 9 images during 
18:44:56—18:56:00 UT, for EUVI/A, 171 A. The observed dimming is represented with hatched areas, while the best-fit model profiles are 
plotted with thick linestyle. 




Fig. 6. Comparison of observed and simulated EUVI base-difference images at 5 times for the observations of STEREO/A 171 A (left 
columns) and STEREO/B 171 A (right two columns). The pre-CME image at 18:36 UT was subtracted in these base-difference images. 
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Table 1. Model parameters 



Parameter 



Best- fit value 



Spacecraft separation angle 
Spacecraft plane angle to ecliptic 
Heliographic longitude A 
Heliographic longitude B 
Heliographic longitude Earth 
Heliographic latitude A,B 
Start time of acceleration 
Start time of modeling 
End time of modeling 
Initial height of CME center 
Initial radius of CME 
Width of CME front 
Acceleration of vertical motion 
Acceleration of radial expansion 
Inital vertical velocity 
Inital expansion velocity 
Maximum density compression 
Corona/CME base density 
Super-hydrostaticity factor 
Mean Temperature (171 filter) 
Temperature width (171 filter) 



a sep = 47.17° 

f3 AB = 2.66° 

I a = -102.4° ±0.9° 

l B = -54.9° ±0.9° 

Ie = -78.8° ±0.9° 

= -8.8° ±0.6° 

to=2008-Mar-25 18:38 UT 

ti=2008-Mar-25 18:36 UT 

ti=2008-Mar-25 18:56 UT 

Ro = 45 Mm 

h = 45 Mm 

wo « 10 Mm 

a h 0.52 ms" 2 

cir — 0.54 m s -2 

Vh ~0m s _1 

vr ^0m s _1 

Qn,max ~ 2.0 

n 6.5 x 10 8 cm" 3 
qx = 1.45 
T 0.96 MK 
AT = ±0.37 MK 



4.6 Comparison of Numerical Simulations with Obser- 
vations 

After we constrained all 14 free parameters (listed in Table 1) 
of our analytical 4D model by fitting some observables, such 
as measured coordinates (Fig. 4) and cross-sectional horizon- 
tal and vertical brightness profiles (Fig. 5), we are now in the 
position to inter-compare the numerically simulated images 
with the observed images, as shown in Fig. 6 for 5 selected 
times, for both the STEREO/A and B spacecraft. The com- 
parison exhibits a good match for the extent of the dimming 
region the the bright lateral rims, both extending over about 
1.5 thermal scale heights above the solar surface. The base- 
difference images of EUVI/A reveal a fairly symmetric CME 
(as the model is by design), surrounded by spherical bright 
rims at the northern and southern CME boundaries (as the 
model is able to reproduce it). 

The model, however, is less in agreement with the ob- 
served EUVI/B images. The extent of the EUV dimming 
region matches relatively exactly, although the observed dim- 
ming region is somewhat cluttered with bright postflare loops 
that appear in the aftermath of the CME, which are mostly 
hidden in the EUVI/A observations. The biggest discrep- 
ancy between the model and the EUVI/B observations is the 
location of the brightest rim of the CME boundary. The 
combination of projection effects and gravitational stratifi- 
cation predict a brighter rim on the west side, where we look 
through a longer and denser column depth tangentially to the 
CME bubble, which is not apparent in the observations of 



EUVI/B. Instead, there is more bright emission just above the 
eastern limb that cannot be reproduced by the model. Appar- 
ently there exists stronger density compression on the eastern 
side of the CME bubble than the model predicts. Another 
inconsistency is the bright loop seen in EUVI/B at 18:51 
UT, which does not match the surface of the modeled CME 
sphere as constrained by EUVI/A. Apparently, there are sub- 
stantial deviations from a spherically symmetric CME bub- 
ble model that are visible in EUVI/B but not in EUVI/A. 
Perhaps a flux rope model could better fit the observations 
than a spherical shell model. These discrepancies between 
the observations and our simple first-cut model provide spe- 
cific constraints for a more complex model (with more free 
parameters) that includes inhomogeneities in the density dis- 
tribution of the CME. 

4.7 Estimate of the CME Mass 

Our model allows us, in principle, to estimate the CME mass 
by integrating the density n e (x,y, z,t) over the entire CME 
sphere, which is of course growing with time, but expected to 
converge to a maximum value once the CME expands far out 
into the heliosphere. A simple lower limit can analytically be 
obtained by integrating the density in the cylindrical volume 
above the footpoint area, 

M C me = rn p J n e (x,y,z)dV C ME ~ m p 7rRlno\ T q\ .(22) 

From our best-fit values Ro = 45 Mm, q\ = 1.45, nO = 
6.5 x 10 8 cm -3 and the thermal scale height of X T = 47 
Mm, we obtain a lower limit of Mcme > 0.47 x 10 15 g. 
However, this CME appears to expand in a cone-like fashion 
in the lowest density scale height, so the total volume and 
mass is likely to be about a factor of « 2 higher. Moreover, 
the mass detected in 171 A amounts only to about a third of 
the total CME mass (Aschwanden et al. 2009b), so a more 
realistic estimate of the total CME mass is about a factor 6 
higher than our lower limit, i.e., Mcme ~ 3 x 10 15 g, which 
brings it into the ballpark of previous CME mass determina- 
tions of this particular event, i.e., tucme = 2.9 x 10 15 g 
from STEREO/COR-2 white-light observations (Colaninno 
and Vourlidas 2009), or m c me = (4.3 ± 1.4) x 10 15 g from 
STEREO/EUVI observations (Aschwanden et al. 2009b). 



5 Conclusions 

We developed an analytical 4D model that simulates the 
CME expansion and EUV dimming in form of a time se- 
quence of EUV images that can directly be fitted to stereo- 
scopic observations of STEREO/EUVI. The dynamic evolu- 
tion of the CME is characterized by a self- similar adiabatic 
expansion of a spherical CME shell, containing a bright front 
with density compression and a cavity with density rarefac- 
tion, satisfying mass conservation of the CME and ambient 
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corona. We forward-fitted this model to STEREO/EUVI ob- 
servations of the CME event of 2008-Mar-25 and obtained 
the following results: 

1 . The model is able to track the true 3D motion and ex- 
pansion of the CME by stereoscopic triangulation and 
yields an acceleration of a w 0.54 m s -2 for both the 
vertical centroid motion and radial expansion during the 
first half hour after CME launch. 

2. Fitting the EUV dimming region of the model to the 
data mostly constrains the coronal base density in the 
CME region (n = 6.5 x 10 8 cm -3 ) and the density 
scale height, which was found to be super-hydrostatic 
by a factor of q\ = 1.45. 

3. The average CME expansion speed during the first 10 
minutes is approximately < v >^ 400 km s _1 , simi- 
lar to the propagation speeds measured for EIT waves 
in the initial phase (i.e., v « 460 km s _1 , Veronig et al. 
2008; v « 475 km s~\ Long et al. 2008; v « 250 km 
s _1 , Patsourakos et al. 2009), and thus the CME expan- 
sion speed seems to be closely related to the associated 
propagation kinematics of EIT waves. 

4. The derived base density, scale height, and footpoint 
area constrain the CME mass, but accurate estimates re- 
quire a more complete temperature coverage with other 
EUV filters (e.g., see Aschwanden et al. 2009b). 

5. The width and density compression factor of the CME 
front are also constrained by our model, but accurate 
values require a perfectly homogeneous CME shell. 

6. While the spherical shell geometry reproduces the 
EUVI/A observations well, significant deviations are 
noted in the EUVI/B observations, indicating substan- 
tial inhomogeneities in the CME shell, possibly requir- 
ing a hybrid model of bubble and flux rope geometries. 

The most important conclusion of this modeling study is 
that EUV dimming can be understood in a quantitative man- 
ner and that it provides a direct measurement of the coronal 
mass loss of a CME. This exercise has shown that the spher- 
ical shell geometry can reproduce a number of observed fea- 
tures of an evolving CME, which constrains the physical and 
kinematic parameters of the initial phase of the CME launch, 
but reveals also significant deviations that require a modifi- 
cation of the idealized homogeneous spherical shell model. 
The method of analytical 4D models with forward- fitting to 
stereoscopic EUV images appears to be a promising tool to 
investigate quantitatively the kinematics of CMEs. Com- 
bining with simultaneous magnetic or MHD modeling may 
further constrain the physical parameters and ultimately pro- 
vide the capability to discriminate between different theoret- 
ical CME models. Full 4D modeling of the initial CME ex- 
pansion may also provide a self-consistent treatment of EIT 



waves and CME expansion (e.g., Chen et al. 2005), for which 
we find similar kinematic parameters. 
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